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Abstract 

We study the pricing and the hedging of claim ijj which depends on the default 
times of two firms A and B. In fact, we assume that, in the market, we can not buy 
or sell any defaultable bond of the firm B but we can only trade defaultable bond of 
the firm A. Our aim is then to find the best price and hedging of tp using only bond of 
the firm A. Hence, we solve this problem in two cases: firstly in a Markov framework 
using indifference price and solving a system of Hamilton- Jacobi-Bellman equations, 
secondly, in a more general framework, using the mean variance hedging approach and 
solving backward stochastic differential equations (BSDE). 
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Introduction 

Models for pricing and hedging defaultable claim have generated a large debates between 
academics and practitioners during the last subprime crisis. The challenge is to model the 
expected losses of derivatives portfolio by taking into account the counterparties defaults. 
Indeed, they have been affected by the crisis and their agreement on the derivatives contracts 
can potentially vanish. In the literature, models for pricing defaultable securities have been 
initiated by Merton [27]. His approach consists of explicitly linking the default risk of a firm 
to its value. This model is a good issue to understand the default risk. However, it is less 
useful in practical applications since it is too difficult to capture the dynamics of the firm's 
value which depends on many macroeconomics factors. In response of these difficulties, 
Duffie and Singleton [9] introduced the reduced form modeling which has been followed 
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by Madan and Unal [26], Jeanblanc and Rutkowski [17] and others. In this approach, the 
main tool is the " default intensity process" , which describes in short terms the instantaneous 
probability of default. This process combined with the recovery rate of the firm represents 
the main tools necessary to manage the default risk. However, we should manage the default 
risk considering the financial market as a network where every default can affect another 
one and the propagation spreads as far as the connections exist. In the literature, to deal 
with this correlation risk, the most popular approach is the copula. This approach consists 
to define the joint distribution of the firms on the financial network with respect to the 
marginal distribution of each firm. In static framework ^, Li [25] was the first to develop this 
approach to model the joint distribution of the default times. But, all computations are done 
without considering the evolution of the survey probability given the available information. 
Thus, we can not describe the dynamics of the derivatives portfolio in this framework. 
In response of these limits on the static copula approach, El Karoui, Jeanblanc and Ying 
developped a conditional density approach [11]. An important point, in this framework, is 
that given this density, we can compute explicitly the default intensity processes of firms in 
the financial market considered. We will follow this approach and work without losing any 
generality in the explicit case where the financial network is defined only with two firms 
denoted by A and B. The intensity process jumps when any default occurs. This jump 
impacts the default of the firm and makes some correlation between them. We assume that 
we can not buy or sell any defaultable bond of the firm B but we can trade a defaultable 
bond of the firm A. We will consider two different cases for pricing and hedging a general 
defaultable claims tjr. the indifference pricing in Markov framework and the Mean- Variance 
hedging approach for the general case. 

In the first case, we will work in a Markov framework. Our aim will be to find, using 
the correlation between the two firms, the indifference price of any contingent claim given 
the risk aversion. This risk aversion will be defined by an exponential utility function. 
We will express the indifference pricing as an optimization problem (see El Karoui and 
Rouge [12]) and we will use the Kramkov and Schachermayer [21] dual approach. Then 
solving this dual problem, we will find the solution of the indifference price. Moreover, 
the characterization of the optimal probability for the dual optimization problem will be 
solved by Hamilton-Jacobi-Bellman (HJB) equations since the defaultable bond price will 
be assumed to be a Markov process in this framework. We will also find an explicit formula 
for the optimal strategy. 

In a second case, we will be interested in a hedging problem using the Mean- Variance 
approach. We will assume that we work in a general setting (not necessarily Markov), then 
we will not be able to others use the HJB equations to characterize the corresponding value 
function. Hence, we will adopt the Mean Variance approach which has been introduced 
by Schweizer in [29] and generalized by many ([30], [13], [22], [8], [1], [23], [14]). Most of 
theses papers use martingales techniques and an important quantity, in this context, is the 
Variance Optimal Martingale Measure (VOM). The VOM, P, is the solution of the dual 
problem of minimizing the L^-norm of the density dQ/dP, over all (signed) local martingale 
measures Q for the defaultable bond price of the firm A. If we consider the case of no jump 
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dues to default, then the bond price process of the firm A is continuous. In this case, Delbaen 
and Schachermayer in [7] proved the existence of an equivalent VOM P with respect to P. 
Moreover the price of any contingent claim tp is given by ¥F{i/;). In Laurent and Pham 
[22], they found an explicit characterization of the variance optimal martingale measure in 
terms of the value function of a suitable stochastic control problem. In the discontinuous 
case, when the so-called Mean- Variance Trade-off process (MVT) is deterministic, Aral [1] 
proved the same results. Since we will work in discontinuous case and since in our case the 
Mean variance trade-off process is not deterministic (due to the stochastic default intensity 
process), we will not be able to apply the standards results. Hence our work will be first to 
characterize the value process of the Mean- Variance problem and then to make some links 
with the existence and the characterization of the VOM. However, we will not really need 
to prove and assume this existence to solve our problem. Indeed, we will solve a system of 
quadratic Backward Stochastic Differential Equations (BSDE) and we will characterize the 
solution of the problem using BSDE's solutions. The main contribution in this part will 
be the explicit characterization of the BSDE's solutions without using the existence of the 
VOM. We will obtain an explicit representation of each coefficients of quadratic backward 
stochastic differential equations with respect to the parameters asset of our model. In 
particular, the main BSDE coefficient will follow a quadratic growth and its solution will 
be found in a constrained space. In a particular discontinuous filtration framework (where 
the asset parameters do not depend on the filtration generated by the jump), Lim [24] have 
reduced this constrained quadratic BSDE with jumps to a constrained quadratic BSDE 
without jump and solved the corresponding BSDE. In the discontinuous filtration due to 
defaults events, we will can not do the same assumption since the intensity processes depend 
on the jumps (the default events). Hence, using Kharroubi and Lim [18] technics, we will 
split the BSDE's with jumps into many continuous BSDEs with quadratic growth and we 
will conclude the existence of the solution using the standard results of Kobylanski [19]. 

Hence, the paper is structured as follow, in a first section, we will give some notations 
and present our model with some results relative to credit risk modeling. Then, in a second 
part, we will study the case of pricing and hedging defaultable contingent claim in a Markov 
framework using indifference pricing. Then in the last section, we will study the pricing and 
hedging problems in a more general framework (not Markov) using mean variance hedging 
approach and solving a system of quadratic BSDEs. 

1 The defaultable model 

In the sequel, we will work in the same model construction as in Bielecki and al. in 
[2] chapter 4. Let T > be a fixed maturity time and denote by (0,F := {J't)[o,T]^^) ™ 
underlying probability space. The filtration F is generated by a one dimensional Brownian 
motion W. Let and be the two default times of two firms A and B. Let define, for 
all t G [0,r]: 

H^ = hr^<t} and Ff = l{,B<,}. (1.1) 
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We define now some useful filtrations and definitions: 

g^ = Ft\J'Hf, g^ = TtVHf and Qt = Ft\J ^ Uf , 

where (resp. H^) is the natural filtration generated by (resp. H^). We will denote 

by G := (^i),,[o,T]. ■= (^/')ie[o,T] := {Q?),^^,^^y 

Definition 1.1 (Initial time). Let -q he a positive finite measure on M^. The random 
times and are called initial times if, for each t G [0, T], their joint conditional law 
given J-t is absolutely continuous with respect to r]. Therefore, there exists a positive family 
idtiu)) te[o T] of ¥ -martingales such that 

Gi(e^^^) = P(T^>e^T^ >0^|J-t)= / / gtiyi,y2)r]{dyi,dy2), (1.2) 

for each 9"^, 9^ € M+ and t € [0,r]. 

Regarding this definition, we make the following assumptions: 
Assumption 1.1. ( Properties of the default times) 

- Processes and have no common jump: P (ja = tb) = 0. 

- The default times ta and Tf^ are initial times. 

Hence, point 2. of the previous Assumption implies that the default times of firm A and B 
are correlated regarding our joint probability density gt (appearing in (1-2)). We now give 
a representation Theorem of our defaultable model. 

Theorem 1.1. (Representation Theorem) Under Assumption 1.1, fori € {A,B}, there 
exists a positive G-adpated process A*, called the ^-intensity of , such that the process 
M* defined by 

Mi = Hl- ! Xids, 
Jo 

is a G-martingale. Moreover, any local martingale ( = {Ct)t>o (admits the following decom- 
position: F-a.s, 

Ct = Co+ f ZsdWs + r UfdMf + /* f/f dMf , Vi > 0, (1.3) 
Jo Jo Jo 

where Z, and are G-predictable processes and W is the martingale part of the 
G-semimartingale W in the enlarged filtration (see [15] for more details about the pro- 
gressive enlargement of filtration and the characterization of the decomposition of any F- 
semimartingale in the enlarged filtration G). 

Proof. The processes A"^ and A^ are given explicitly since we assume that and are 
initial times and given our conditional law G. Moreover in Proposition 1.29, p54 in [31], the 
author follows the proof of the representation Theorem of Kusuoka (representation theorem 
when the default times are independent of the filtration F) to construct the proof when 
default times are initial. □ 
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1.1 Dynamic of the Bond 

In our model, the traded asset will be the defaultable bond of the firm A. Using the 
decomposition (1.3), we represent the dynamics of this defaultable bond in the enlarged 
filtration G as in Corollary 5.3.2 of [2]: 

= fitdt + afdUt"^ + afdMf + atdWt, (1.4) 

t- 

where fi,a^,a^ and a are G-predictable bounded processes. Therefore, given an initial 
wealth X > 0, if we assume that investors follow an admissible strategies vr, which is 
represented by a set A of predictable processes vr such that 



E 



T 

' < +CXD. (1.5) 



Then we can define the dynamics of the wealth process, started with an initial wealth x at 
time t = and following a strategy vr, X^''^ based on the trading asset D"^ by 

dXf = ^t— ^ = vrt [iitdt + cj^dM^ + afdM^ + atdWt] . (1.6) 

Note that since all the coefficients in the dynamics of the wealth process are bounded, then 
for any tt € .A we have that (1.5) implies: 



E 



. , t I 

0<t<T 



< +00. 



1.2 The Defaultable claim 

We now introduce the concept of defaultable claim and give some explicit examples. 

Definition 1.2. A generic defaultable claim ^ with maturity T > on two firms A and B 
is defined as a vector 

i;:= (X^,X^,Z^,Z^,T^,T^) 

with maturity T such that: 

- r', i € {A,B} is the default time specifying the random time of default of the firm i 
and thus also the default events {r* < t} for every t € [0,r]. It is always assumed 
that T* is strictly positive with probability 1. 

- X"^ is the promised payoff which represents the random payoff received by the owner 
of the claim ip at time T, if there was no default of firm A prior to or at time T. 

- X^ is the promised payoff which represents the random payoff received by the owner 
of the claim tp at time T, if there was no default of firm B prior to or at time T. 

- Z'*, i E {A^B}is the recovery process which specifies the recovery payoff Z^i received 
by the owner of a claim at time of default of the firm i, provided that the default 
occurs prior to or at maturity date T. 
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We can introduce now the payoff at time T of this defaultahle claim, which represents all 
cash flows associated with {X^^, , Z^, ,t^,t^). We will use also the notation for 
this payoff. Formally, the payoff process ip is defined through the formula by 

= X\,A^T} + X''l{,B^T}+ f ZtdH^+ r ZfdHf. (1.7) 

JO JO 

As an example, we can have a defaultable claim which only gives a terminal payoff of 
if no default occurs before time T. Hence, we will not receive money if one of the firms 
makes default. So our defaultable claim is given by 

We can also have a defaultable claim which gives an amount of money with respect to the 
default time of the firm B and gives a recovery amount if the firm A makes default 

1 2 r A 

2 Hedging defaultable claim in Markov framework 

Let consider ip G Qt a bounded defaultable claim as defined in Definition 1.2, which 
depends on the default times of the firm A and of the firm B. Our aim is to find the 
best hedging and pricing of ^p with respect to these defaults times. 

Assumption 2.2. We assume that ji^a^ and the intensity processes \"^,\^ are 
deterministic bounded functions of time, and . 

Remark 2.1. Under Assumption 2.2, we have that {D^,H^,H^) is a Markov process. 

We assume that the risk aversion of the investors is given by an exponential utility function 
U with parameter 6, given by 

U{x) = — exp(— (5x). 

Therefore, to define the indifference price or the hedging of ■0, we should solve the following 
equation: 

u^{x + p) = u'^{x), 
where functions and are defined by: 

u^{x) = supE [-exp(-5(X^'^ - -0))] and u°{x) = supE [- exp(-(5X^'^)] . (2.8) 

2.1 The dual optimization formulation 

To deal with the problems (2.8), we use the duality theory developed by Kramkov and 
Schachermayer in [21]. In fact this theory allows us to find the optimal wealth at the 
horizon time T and the optimal risk- neutral probability Q*. Let recall now some results 
about the dual theory. 
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Theorem 2.2. [Kranikov and Schachermayer, Theorem 2.1 of [21]] 

Let U be an utility function which satisfies the standards assumptions and consider the 
optimization problem: u{x) = supj^g_4E [[/(X^'^^)], then the dual function of u defined by: 



v{y) = sup{n(x) — xy], u{x) = \iii{v{y) + yx} 



is given by 



viy) = inf E ( y 



y 



dF 



(2.9) 



where V represents the dual function of U and M.^ represents the set of all risk-neutral 
probability measures. 

Moreover, there exists an optimal martingale measure Q* which solves the dual problem and 
we have that the optimal wealth at time T is given by: 



, where v is defined s.t. E^ 



The function I represents the inverse function of U' and represents the Radon Nikodym 
density of Q* with respect to P on Qt- 

We can apply this result to solve our optimization problem (2.8). We will resolve only the 
case ip ^ 0. Indeed the particular case 1^ = could be obtained as a particular case of these 
results. We obtain an analogous result of Delbaen and al. Theorem 2 in [5], given by the 
following proposition: 

Proposition 2.1. Let Q* be the optimal risk-neutral probability which solves the dual prob- 
lem 



inf 



H{q\F) - 6E'^{t/j) 



(2.10) 



then the optimal strategy vr* € ^ solution of the optimization problem (2.8) satisfies: 



1 



1 



-- In (^Z^ j + = X + - In + TTldDf, 

where i7(Q|P) represents the entropy of Q with respect to P (j,.e. ES 
is a non negative constant. 



log 



(2.11) 

and y 



Proof. The proof is based on the Theorem 2.2. Firstly, to match with assumptions of this 
theorem in the case ^ ^ 0, we change the historical probability. Let define 

Gt E[exp(5'0)] Treyl 



then setting c = E [exp((5'0)], we get 
n^(x) 

1 



supE [-exp(-(^(X^'" - V'))] = supE^"* [-cexp(-5X^'")] 



supE'' 



exp 



log(c) + X. 



X,TT 

T 



= sup E^"^ 


exp 1^ 


tt€A 





log(c),7r 

'■T 
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Hence by the definition of u^{x), we obtain tliat u'*' \ x 



■7^ 



I ln(c)) = u^{x). Tlien using the 
Theorem 2.2, the dual function of is given, for all y > 0, by : 



^(y) = inf E 



Vly 



(2.12) 



where 



y 

V{y) = sup{[/(x) — xy} = sup{— exp(— fe) — xy} = — 

x>0 x>0 



1 



Using this expression into (2.12) gives, after straightforward calculation, an explicit expres- 
sion of the dual function which is given by 

^(y) = V{y) + ^ ln(c) + I inf \h{Q\F) - 5E^(V')" . 
Q^M" I 

Since Q* is the optimal risk-neutral probability which is solution of (2.10), we deduce that 
the optimal wealth at time T of the optimization problem (2.8) is given by 



where y is defined such that E*- 



X — i ln(c) and I is equal to —V'. 



Moreover from Owen in [28] , we can deduce that there exists an optimal strategy ir* & A 
such that: 



z: 



y- 



T 



z: 



T J 



x-^ln(c) + ^ T^ldDt 



In our case, since we work under the case of an exponential utility function with parameter 
5, we have that 



We finally get 



1 



X - - ln(c) + I TT^dD. 



T 



S JO 

which concludes the proof of this proposition. 



i>n(|)-i,o.u: 



T I + '4' - ^ln(c). 



□ 



2.2 Value function of the dual problem 

We are now interested in solving the dual problem. Firstly, let us consider the same 
problem with a different set of probability measures like = Q, where Q represents the 
set of all probability measures Q ^ P. Then the value function is given by the entropy of ip 
with a parameter 5. But since we work in a more restricted set of probability measures 
which represents the set of all risk-neutral probabilities, the value function is then more 
difficult to precise. Indeed, to characterize the value function, we begin by defining the 
set M^. Hence, let Q G and define Z'^ to be the Radon Nikodym density of Q with 



respect to P. Considering the non negative martingale process Z^ 



E 



Z^lQt 



and using 
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representation Theorem 1.1 imply that there exists predictable processes p"^ and which 
take their values in C = (—1, +00) and a predictable process p which takes its values in M 
such that for all t G [0, T], 

dZ^ = Zl [pUm^ + p^Mf + ptdW,) . 

Since Q is in it is a risk-neutral probability, then ZD^ is a local martingale. This 
implies by Ito's calculus the following equation: 

pt + P^o^X^ + pf af Af + ptat = 0. (2.13) 
Remark 2.2. We notice that the process p depends explicitly on the values of p^ and p^ . 
Therefore using equation (2.13), (2.10) can be view as find p^ and p^ which minimize: 



inf ¥1^ 



\n{Z^)-5i} . (2.14) 



This is the dual problem we would like to solve. We make now an assumption on the 
decomposition form of our defaultable claim ip. 

Assumption 2.3. The defaultable claim ip ^ Qt is given by 
where g and f are two bounded continuous functions. 

Remark 2.3. 1. We choose to take a defaultable claim which depends only on the de- 
fault time of the firm B. However, we could have been take a defaultable claim which 
depends on the default time of the firm A too. The calculus would have been longer 
hut the results would have been the same. 

2. Moreover, taking a defaultable claim depending only on the default time of the firm B 
has an economic sense. Indeed, our traded asset is the defaultable bond of the firm A, 
so it is justified to take payoff g and f function of D^. Therefore if we see the firm 
B as an insurance company which covers the firm A, then the default of B means the 
counterparty default risk. 

Proposition 2.2. Under Assumption 2.3, the value function of the dual problem (2.14) is 
given by: 

V{t,DtH^,H^) := inf M T j{s, pt p"^ ,Dj)ds - 5g{D^)l<,B^T} Dt,H^,H^ 



^ j{s, pf,pf, Df)ds - 5g{D^)l{,B^T} 



(2.15) 

where the function j is defined by: 

j{s,pf,pf,Dt)= K[{i + pl)Hi + pl)-pl]-s{i + pf)xff{Df) + ^-pl 

ie{A,B} 

(2.16) 
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Proof. The proof is based on the Ito's formula. The dynamics of ln(Z'^) under Q is given 
by 

dln(Z?)= Yl pldMi+[ln{l + pi)-pl\dHi + ptdWt-^pUt. 

i€{A,B} 

Using Girsanov theorem, the processes defined for all i S {A, B} by 



Ml = Ml- / p\\\ds and Wt = Wt 



Psds 





are Q-martingales. Hence, we obtain that 

ln{Z^) -5^=rY. + P^) 1^(1 + (^t) - P\]dt -S \ r f{Df^ )dHf + g{D^){l - +/ 

■'%&{A,B} LJo J Jo 

where is a Q-martingale. Then, we can rewrite the dual problem using the last expres- 
sion: 

r /■'^ 

ln(Z^)-5V = inf 



-pidt+M^^ 



inf E*^ 



j{s,pf,pf,Df)ds-6{l-H^)giD^ 







where j is given in (2.16). Since by Remark 2.1, the process {D"^, H^, H^) is a Markov 
process, then using the standards results of [4], the value function of the dual optimization 
problem is given by: 



^ jis, pf,pf, Df)ds - 5g{D^)l{,B^T} 



□ 



Proposition 2.3. Let z = {x,h'^,h^) and h = {h^,h^), then the value function of the 
dual optimization problem is solution of the following Hamilton- J acohi- Bellman equation: 

^it,z) + l^{t,z)a\t,z) +Jnl^£^A^pBVit,z)+j{t,pt,pf)} = 0, V{T,z) = g{x){l - /i^) 



(2.17) 



where 



CpApBV{t,z)= 

i€{A,B} 



dV 

dz 



{t,z)a\t,z) + {V{t,z')-V{t,z)) 



{l + pDX\t,h), 



and = [x{l + cr^it, z)), + a*, + \ — a*) where = 1 and = 0. Moreover given 
the value function, the optimal strategy satisfies: 



Pt 



Df.a{t,z) 



where the process p is explicitly given with the optimal control p^ , i G {A, B}, see the relation 
(2.13). 
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Proof. From Proposition 2.2, we find that the value function of the dual optimization prob- 
lem is given by (2.15). Since [D^, H^, H^^ is Markovian under P and that the risk neutral 
probability measure Q depends on the control {p^,p^), we can apply the same method as 
in [4] section 3.2 and 3.3. So, using now Hamilton-Jacobi-Bellman (HJB) equation we get: 

Vit, Df, H^, ) = inf [ jis, pt pf,D^)ds + Vit + h, H^_,^, H^+hW, H, 
Then the value function solve the HJB equation (2.17). 

We are now interesting in finding the optimal strategy given the value function. Let recall 
that from Theorem 2.2, the optimal risk- neutral probability and the value function exist. 
Let define p^,p^ and p be the optimal density parameters. Since p^ and p^ are optimal for 
the HJB equation, assuming cj(t, z) 7^ , using first order condition, we find for i G {^4, B}: 



B 



dV crHf z) 

{V{t, z^)-V{t, z))-xa\t, z) — {t, z)+ln(l + p]) ^pt 

ox cr[t,z) 



X\t, h)=6{l - a')f{x)\\t, /i). (2.18) 



Then using the HJB equation (2.17) and the relation (2.18), we obtain the following relation: 

P1 



\pI+ plX%h) 

ie{A,B} 



(1 + Pi)^Fr^~Pt + n ^)^V2(t, z) + — (t, z). (2.19) 



iG{A,B} 



a{t,z) 



2 dx"^ 



dt 



Let recall the Ito's decomposition of the process ln(Z'^*): 



\n{Z^')= r[ptdWt + \pldt] + r H^ + mHl--^^\\t,h)\. 

Then using equations (2.18) and (2.19), we get an useful and more explicit decomposition 
of the process ln(Z^ ): 

-{t,zt){D^-) a^{t,zt)dt 



ln{Zl 



T 



2 9x2 







{t,zt)dt+ / ptdWt 



dt 



E 

i&{A,B} 
T 



{V{t, z\) - y (t, zt)) - Df. a\t, zt)^{t, z) 







dHl 



+ 



E ^-^:^Pt[dHl-{l + m{t,ht)]+ r5f{Df^)dH^, 



ie{A,B} 



a{t,zt) 



where zt = (Df, H^, H{^) and ht = (Hf, H^). The Ito's decomposition of V{T, D^, H^, 



gives: 



ln(Z^ 



T 



Pt 



a{t,zt)dWt+ Y ^\t.zt)dMi 

i(i{A,B} 

dx 



a{t,zt) 

ViT,D^,H^,H^) + V{0,D^,HlH^)+ [ 

Jo 



{t,zt)dDf. 



Since 

V{T,D^,H4,H^) = -6g{D^){l-H^) 
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and 



we get: 



ln(Zf ) - <5V = ^(0, D^, H^, 



T 



Pt 



dV 



Finally, from the definition of the value function, we have 



ln(4 



and using the fact that E*^ 
Theorem 2.2), we deduce that 



^ ln(c) where X^'' 



dDi 



-lln ( 



(see 



Hence, we conclude 



]ln(Zf 



+ V-Jln(c)-iln(| 



X — - ln(c) 





V{0,D^,H^,H^ 



and we find 



-^ln(Z^' 



+ ^ = x + ^In f^) + 

\d/ ,IQ 



-6x — In ( - 



Pt dV 



Therefore from equation (2.11), we obtain the expected result. 



dDi 



□ 



In conclusion, we found that since we can characterize the optimal probability for the 
dual optimization problem using Kramkov and Schachermayer Theorem, we can charac- 
terize the HJB equation solution of our dual problem. This allows us to find the optimal 
strategy for the primal solution for a defaultable contingent claim ip. Therefore we can 
find for -0 = and -0^0 the optimal strategy and deduce the indifference price p of a 
defaultable contingent claim solving the equation u^{x +p) = u^{x). 



3 Generalization of the hedging in a general framework: Mean- 
Variance approach 

In this part, we assume that we work in a more general setting (not necessarily Markov), 
then we cannot use the HJB equation to characterize the corresponding value function. To 
solve our problem, we will use the Mean Variance approach. It is a well-known methodology, 
introduced by Schweizer in [29], to manage hedging in general case. An important quantity 
in this context is the Variance Optimal Martingale Measure (VOM). The VOM, P, is the 
solution of the dual problem of minimizing the L^-norm of the density over all (signed) 
local martingale measure Q for D^. Let recall now the Mean- Variance problem: 

V{x) = minE [(X^'^ - V)^] • (3.20) 
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If we assume that G = F (in this case we do not consider jump of default), then the process 
is continuous. In this case Delbean and Schachermayer in [7] proved the existence of 
an equivalent VOM ¥ with respect to P and the fact that the price of ip is given by ¥F{iIj). 
In the discontinuous case, when the so-called Mean- Variance Trade-off process (MVT) (see 
[29] for definition) is deterministic, Aral [1] proved the same results. Since we work in 
discontinuous case and since the Mean Variance Trade-off process is not more deterministic 
(due to the stochastic default intensity process), we cannot apply the standard results. 

Remark 3.4. Indeed, in this part we do not assume anymore that intensity processes 
and to be deterministic. We take general stochastic default intensity processes. But we 
assume that default times and are ordered, < and that the (Ti)- hypothesis 
holds. A financial interpretation of this assumption could be the counterparty risk. Indeed, 
the firm A could be a bank (counterparty) and the firm B its company insurance which 
covers its default. 

So our work is firstly to characterize the value process of the Mean- Variance problem 
using system of BSDE's. Secondly, to make some links with the existence and the char- 
acterization of the VOM in some particular cases. Thirdly, to prove the existence of the 
solution of each BSDE and to give a verification Theorem. 

We begin by recalling some usual spaces: 

• For s < T, S°°[s,T] is the Banach space of M- valued cadlag processes X such that 
there exists a constant C satisfying 



• For s < T, 'H'^[s,T] is the Hilbert space of M-valued predictable processes Z such that 



• BMO is the space of G-adapted matingale such that for any stopping times < a < r < T, 
there exists a non negative constant c > such that: 



then M = Z.W G BMO, to simphfy notation we write Z G BMO. 

Definition 3.3 (i?2(lP) condition). Let Z be a uniformtly integrable martingale with Zq = 1 
and Zt > 0, we say that Z satisfies the reverse Holder condition R2{¥) under ¥ if there 
exists a constant c > such that for every stopping times a, we have: 



^\\s°°[s,T] •= sup \Xt\ < C < +00. 



t€[s,T] 




E [[M]r - [M]^- \g^] < C, 
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3.1 Characterization of the optimal cost via BSDE 

On our problem of mean-variance hedging (MVH) (3.20), the performance of an ad- 
missible trading strategy vr G ^ is measured over the finite horizon T for an initial capital 
x > by 

J^(T,7r) = E[(X^'^ - V)^]. (3.21) 



We use the dynamic programming principle to solve our mean variance hedging problem. 
Let first denote by A{t, u) the set of controls coinciding with v until time t € [0, T] 

A{t, v) = {-KeA: vr.At = v.m}- (3.22) 
We can now define, for all t € [0,T], the dynamic version of (3.21) which is given by 



J'^(t,7r) = essinf E 



(3.23) 



Let recall now the dynamic programming principle given by El Karoui in [10]. 

Theorem 3.3. Let S he the set of G- stopping times. 

1. The family {J^ {t,v),t € 5, G A} is a submartingale system, this implies that for 
any u ^ A, we have for any a < t, the following submartingale property: 



E 



a.s. 



(3.24) 



2. ly* £ A is optimal if and only if {J"^ {t,v*),t S 5} is a martingale system, this means 
that instead of (3.24), we have for any stopping times a < t that: 



E 



a.s. 



3. For any u G A, there exists an adapted RCLL process J'^{v) = {J'^{v)t) 
is right closed submartingale such that: 

= J^{t, I'), IP — a.s, for any stopping time r. 



0<t<T 



which 



We search as in Lim [23] a quadratic decomposition form for jf as 



(3.25) 



such that G is a non-negative G-adapted process and Y, ^ are two G-adapted processes. So, 
we assume the quadratic form (3.25) of the cost conditional J"^ with respect to the wealth 
process and we will use the Theorem 3.3 to characterize the triple (0, Y, ^) as solution 
of three BSDEs. We will verify in the section 3.2 that the assumption of the quadratic 
decomposition form, the optimality and admissibility of the founded optimal strategy are 
satisfied. 
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So, let TT € ^ be an admissible strategy, by representation Theorem 1.1, we have that 
the triplet (0,y, ^) need to satisfies the following BSDEs: 



dQt 

dYt = -g^Yt, Ut^, Uf, Zt)dt + UfdMf + UfdMf + ZtdWt, 
dit = -gli^t, et ef , Rt)dt + efdM^ + ef dM^ + RtdWt, 



Gt = 1 



(3.26) 



with the constraint that 6t > 5 > 0, for some non negative constant 6, for all t € [0, T]. 
The processes are G-predictable. Hence, we can use Ito's formula 

and integration by part for jump processes to find the decomposition of J'^{'k). Let recall 
that for any semimartingale S and L, we have that 

d{StLt) = St-dLt + Lt-dSt + d[S, L]t. 

In our framework since a jump comes from defaults events we get 

d[S,L]t = {S',L')t+ Yl ^SlALldHi. 

i£{A,B} 

Applying these results for 5 = L = (X^'^ — Y) gives: 



diX"^'^ - Yf^ = 2{X^,f - Yt- 



{TTtfit + gt)dt + {T^tal - Ui)dMi + {TTtat - Zt)dWt 

i£{A,B} 



+ {atTTt-Ztfdt+ J2 M-UlfdHl 

iG{A,B} 

Secondly take 5" = 6 and L = (X^''^ - Y f , let define K := (X^-'^ - Y), we find: 

d{eK^)^ = 2Kt-et- i'KtfJ.t + 9t)dt+ Yl {7Ttal-Ui)dMi + {7Ttat-Zt)dWt 
I ie{A,B} 

+ e^- (atTTt - Ztfdt + Y - UifdHl - Qt- Kl gjdt 



i£{A,B} 

+ et-Kl Y &idMi + (3tdWi 

_ie{A,B} 

+ E [M-Ulf + 2K,-{7rtal-Ul 

i£{A,B} 



+ 2Kt-et-{'^tC7t- Zt)Ptdt 

Ol@t-dHl 



Using this decomposition, we can write explicitly the dynamics of J^(vr) for any tt £ A, 
djf {t:) = dM^ + dV^ where M/^ is the martingale part and the finite variation part of 



4>. 
t ■ 



djf'in) = dM^ + e^- [TT^at + 27rt{btKt + q) + 2Kt{g^ - ut) - K^gj + vt] dt - g^dt 



(3.27) 
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where processes are defined respectively by: 

iG{A,B} iG{A,B} 

2 



Ct 



ie{A,B} ie{A,B} 

and ut = f3tZt+ ^ ^^^t^t- 



(3.28) 



ie{A,B} 



Using now Tlieorem 3.3, we have that, for any tt ^ A, the process J'^(vr) is a submartingale 
and that there exists a startegy tt* & A such that J^{'k*) is a martingale. This martingale 
property implies that we should find it* such that the finite variation part of J^'{'k*) van- 
ishes. Since the coefficients g^^g^ and do not depend on the strategy vr, using the first 
order condition, we obtain: 



7r+ 



btKt + Ct 
at 



t < T where Kt = - Yt. 



(3.29) 



Therefore, substituting this explicit expression of the optimal strategy in (3.27), we obtain: 



{btKt + Ct) 



at 



+ 2Kt{gt-ut)-Ktgl + vt 



dt - gfdt 



62 



-Kngl + "-^]+2Ktig^-ut-^-^ 



Then setting gj + ^ = 0, g'^ - ut - ^ = and {vt - ^)0t- - gf = 0, we find that our 
coefficients g^,g^ and g"^ are given by: 



at 



at 



dt +{{vt- -^)et- - gt 
at 



g]{@t.etof,Pt 



{A,B} ^t'^t \\ + ^t^t 



gt{Yt,U^,U^,Zt) 



+ E ^MAj + AZi 

jg{A,B} 



gt{it,etef,Rt) = Qt 



z't+ {uif{i + em 



Zt'yt + T.;e{A,B}4ui{^ + 0DK 

^e{A,B} <^f + E^e{A,B}i^ + 0l){c7lf\l 

Moreover the solution of the optimization problem (3.20) follows the quadratic form: 

V{x) = Go{x-Yof + Co. (3.30) 
We are now interesting in the proof of the existence of the solution of each BSDE. 
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Remark 3.5. (Existence of the third BSDE) 

1. If we find the solution of the first BSDE (0,9^,6^,(3) e S°°[0,T] x S°°[0,T] x 
S°°[0, T] xBMO, with the constraint O > 6 > and the second BSDE {Y, U^,U^,Z) e 
S°°[0,T] X S°°[0,T] X 5°°[0,T] x BMO then the solution of the third is given by: 



6 = IE 



ds 



t < T. 



Then |(^| € 5°° and from representation Theorem 1.1, we deduce that the martingale 
part M of i: 

Mt= f V e\dMi + / RsdWs 
■^0 i^{A,B} 

is BMO. Moreover from Lemma 3.1, ande^ are bounded. Therefore {^,e^,e^ , R) G 
5°°[0,r] X 5°°[0,r] X S°°[0,T] X BMO 

2. In the complete market case, we have that the tracking error ^ = since the hedging 
is perfect. 

Now, we give the Theorem which proves the existence of the solution of the first 
quadratic BSDE. 

Theorem 3.4. There exists a vector [0,9^ ,9^ , fi) G 5°°[0, T] x 5~[0, T] x 5°°[0, T] x BMO 
solution of the quadratic BSDE 

^ = -5t (Gt, 6^ 9f,f3t)dt + 9tdM^ + 9fdM^ + (3tdWt, Ot = 1. 

Moreover there exists a non negative constant 5 > such that 0t > 5 for all t € [0,T]. 
Given (G, 6*^, 9^,p), we can prove the existence of {Y, U^, , Z) G 5°°[0, T] x <S°°[0, T] x 
5°°[0,r] x BMO solution of the second BSDE {g^,t{;) and {i,e^,e^,R) G 5°°[0,r] x 
5°°[0,r] X S°°[Q,T] X BMO solution to the third BSDE (^^,0). Moreover, given this triplet 
solution (0,y, ^) of our system of BSDEs (3.26), the solution of the our optimization prob- 
lem (3.20) is given by: 

y(x) = Go(x-yo)' + eo. 

The proof of this Theorem will be given in the sequel in section 3.4. 



3.2 Verification Theorem 

Given the solution of the triple BSDEs in their respective spaces (Theorem 3.4), we 
have to verify that the assertions defined in Theorem 3.3 hold true (i.e. the submartingale 
and martingale properties of the cost functional J are satisfied and the strategy vr* defined 
in (3.29) is admissible). Moreover, we will prove that the wealth process associated to vr* 
exists (satisfies a stochastic differential equation (SDE)). 

We begin by proving the existence of a solution of the SDE for the wealth process 
associated to vr*. 
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Proposition 3.4. Let vr* be the strategy, given by (3.29), then there exists a solution of 
the following SDE: 

dXf ^* = vr; [^itdt + a^dM^ + afdM^ + atdWt] with = x. (3.31) 

Moreover, vr* is admissible (i.e. tt* ^ A). 

Proof. The proof is divided in three steps. Firstly, we prove the existence of a solution of 
the SDE satisfied by the wealth process associated to vr*. Secondly, we prove the squared 
integrabihty of this wealth at the horizon time T. Finally, we prove the admissibility of the 
strategy vr*. 

The existence of the solution of the SDE for the vi^ealth process: Plotting the ex- 
pression of TT* given by (3.29) in (3.31) gives 

dXf = {htX^'^' + -ct)dt + ((i?Xi"'"* + et)dMf' + {dfX^^^' + ef)dM^ + {dtX^^^' + et)dWt 

(3.32) 

where the bounded processes are given by 

ot = /it, ct = [bt h ct)fit, dt = (Tt 

at at at 

Y — b ■ ■ Y 

et = {bt— + ct)at, 4 = — -(^u 4 = {bt— + ct)crj, 
at at at 

and processes a, b c are defined in (3.28). We recall, now, that the solution of the 
SDE: 

d4>t = cl)t- [btdt + dfdMf + dfdMf + dtdWt] with (l)o = x 
is given explicitely by 




Setting X^'^ := Ltcpt with 

dLt := qtdt + ifduf + ifdUf + kdWt, Lq = 1, 

we find by integration by part formula that dX^'^ = (pt-dLt + Lt-dcpt + d[(j),L]t. 
Hence, 



dXf = [btdt + d^dMf + dfdMf + dtdWt] + qt - Yl 

i(i{A,B} 

+ (l^t-ll{^ + dl)dMi + cl)t-ltdWt + lt(l)t-dtdt. 

ie{A,B} 

Therefore from equation (3.32), we find, for i G {A,B}, that = (/)j-/j(l + dl), 

et = (pt-h and q = 4)t- [qt — X^igjA B} '^i't'^*) • deduce that the process L is 
defined by: 



Lt = l + 



1 



EAi pi 
(1 + 4) 

i€{A,B} ^ 



ds + 



/* ^dWs + f 

Jo 9s- Jo 



ie{A,B} 



[i + dD 
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and Xt 



tLt is a solution of the SDE (3.31). 



Squared integrability of the strategy vr*: Let us prove first that X^''^* € ?^^[0,T] and 
X^'""* G L2(0, Gt). We recah that jfin*) = Gt(Xf - Ytf+^t is a local martingale. 
Therefore, there exists a sequence of localizing times (Ti)-gj^ for jf such that for 
t< s <T 



E 

From Remark 3.5, we have: 

E l^tAT, - Co] 



-E 



tATi / 2 







a 



where u, c and a are defined in Proposition 3.26. Since a > 0, we have: 

rtAT, 



E 







Moreover, since there exists a constant 6 > such that Qt > 5 and the process v is 
non negative, we can apply Fatou lemma and we find when i goes to infinity that 

ft 



5E 



{XI 



< E 



YA 



<Qq{x-YqY + E 



Therefore Z is BMO and the process e\U' € S°°[0,T] for i = {A,B}. We conclude 
V G 'H'^[0,T]. Hence, we have: 

Since Y G 5~[0,r], then we get the expected results: X^''^* G ■^^[Cr] and X^'''* G 

Admissibility of the strategy vr*: Let now prove that the strategy vr* G 'H'^[0,T]. Ap- 
plying Ito's formula to (A^''^*)^ we get diXt'^^'f = 2X^f* dX'^'^ + d[A^'^*]t, then 
there exists a sequence of localizing times (Tj) -gp^ such that for all t < s < T: 



+ E 



TAT' 



< E 



x,7r 
TAT 



. 2 



2E 



T^Ti 



vr^^.Af' ds 



(3.33) 

Setting = o"^ + (cr^)^A^ + (cr'^)^Af (ET"' is the so called mean variance trade-off 
process), since the processes cr, cj*,A* are bounded, there exists a constant K such 
that > -fC. Then, we obtain 

-27r>,Ar* < -||A-''-*|V.l' + f Kl', < s < T. 



Therefore, combining this inequality with (3.33) gives 



TAT* 



wirKids 



< E 



(X 



TATi^ 



-FE 



TAT 



2 * 2 9 



K 

+—E 
2 



TAT 



|7r*pds 
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Applying Fatou's lemma, when i goes to infinity, we get: 







< E 



+E 



T 



I/is I ds 



K 
+— E 
2 



7r„ as 



uo 



Therefore since K'^ > ii', we finally obtain: 



— E 
2 



- rT ^ ■ 




1 W*A ds 


< E 


Jo 





(X 



a;,7r 
T 



+ E 









Since fj, is bounded, X^''^* G 'H'^[0,T] and X^'''* G L2(17,^t), we conclude vr* G 
■H^[0,T], so vr* is admissible. Note that this condition implies that X^''^ G 5^[0,T] 
since all coefficients of the asset are bounded. 

□ 

We now prove the submartingale and the martingale properties of the cost functional. 

Proposition 3.5. For any tt £ A, the process J'^(vr) is a true submartingale and a mar- 
tingale for the strategy vr* given by (3.29). Moreover the strategy vr* is optimal for the 
minimization problem (3.20). 

Proof. Firstly, we prove the submartingale and the martingale property of the cost func- 
tional then secondly we prove that the strategy vr* is optimal. 

First step: Let recall that for any tt £ A, the process J^{t^) is a local submartingale and 
for TT*, J^{^*) is a local martingale. Therefore, there exists a localizing increasing 
sequence of stopping times (Tj)^gf^ for such that for t < s < T: 

Jf^^^ (vr) < E [ J,AT, (vr) I ] and jf^^^ (tt* ) = E [ J.^t, (vr* ) | ] for any vr G ^. 

(3.34) 

Moreover, for any vr G jf {t^) = ^ti^t'^ ~ ^) +'Ct where , Y and ^ are uniformly 
bounded and X^^'^ G 5^[0, T]. Hence, taking the limit in (3.34) when i goes to infinity 
and applying dominated convergence Theorem, allow us to conclude. 

Second step: For any vr G w4, we have from the submartingale property of J'^(vr) and the 
martingale property of J^{'k*) that : 



E 



< Jo^(^) = eo(x-yo)' + eo = E 4'(vr*) 



Finally, vr* is the optimal strategy for the minimization problem (3.20). 



□ 



3.3 Characterization of the VOM using BSDEs 

Theorem 3.4 leads us to construct the VOM in some complete and incomplete markets. 
We can find also the price of the defaultable contingent claim ijj via the VOM. We consider 
three different cases: 
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i. Complete market (where we assume G = F and G = M^). 

ii. Incomplete market (where we consider only the case G = F VM^). 

iii. Incomplete market (where we consider the case 

Remark 3.6. - The case iii. corresponds to the more general case where the model 
depends on the market information (i.e. the filtration ¥) and the defaults informations 
of the firms A and B. Indeed, in this set up, the model depends on the default time of 
both firms. An economic interpretation of this case is a market with two firms where 
A is the main firm and B its insurance company. Then, the main firm A could make 
default and cause a default of its insurance. Then it is what we call a counterparty 
risk. 

- The case ii. corresponds to a particular case where our model depends only on the 
default time of the firm A. In fact, in this set up, the model depends on the market's 
information and the possible default of the firm A. It can be view as a particular case 
of iii. with condition = oo (i.e. no possible default of firm B). 

- The first case in i., if G = ¥, corresponds to a model which depends only on the 
market information and not to the possible default of firms A and B. In a economic 
point of view, it is a simple model without default. In the second case, G = M'^, 
the model depends only on the possible default of the firm A and non more on the 
information given by the market. 

Remark 3.7. We have explicit solution of the VOM with respect to the process Q in the 
first two cases. 

3.3.1 Complete market 

If we assume that G = F (we do not consider the default impact of firms A and B on 
the asset's dynamics of the firm A) or G = (we do not consider the market noise) then 
our financial market is complete. Hence, the VOM is the unique risk-neutral probability 
and its dynamics can be found explicitly. Our goal in this part is then to find the solution 
of the triple BSDEs given the VOM P. 

Proposition 3.6. Let P be the VOM (the unique risk-neutral probability) and let Zt be 
the Radon Nikodym density of P with respect to P on Qt- We denote Zt = ¥, ^Zt \ Qt] , then 
for all t <T, we have that 

Of — 1 n , T • 

E [Z^lGt] 

Moreover, for all t £ [0, T], we have that Yj = E and S^t = 0. 

Proof. We will consider the two cases G = F and G = H"^. 

First case: Let consider the case where G is equal to F and let the process L be defined 
by the stochastic differential equation given by dLt = L^-ptdWi where pW € BMO. 
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Using Ito's formula we obtain: 



d (^) = §1 [(2P* - f3t)dWt + (/32 + gl - 2f3tpt + pi)dt\ 



{2pt - (3t)dWt +[{fit- Ptf -{- + Pt) 



dt 



{2pt - f3t)dWt + ( i-pt - ^){2/3t + pt + ^)\dt 



Then, if we set, for all t <T, pt := and used the bound conditions of a) 
and the BMO property of /3, we obtain that the process is a true martingale. 
Therefore we get: 



E 



t < T. 



Since Qt = 1, we find the expected result. Moreover, we obtain that L = Z which is 
the Radon-Nikodym of the unique risk- neutral probability and = —^Zt , gf = 0. 
Finally, Yt = E [iplGt] and = 0, t < T. 

Second case: Let now consider the case where G is equal to HI and let the process L be 
defined by the stochastic differential equation given by 

dLt = Lt-pfdMt'^, 



where p^M^ S BMO. Using Ito's formula we find: 





'({I + Pff 


&t- 


A 1 + ^*^ 




'({i + pff 


Ot- 


A 1 + ^*^ 




'({i + pff 




A 1 + ^*^ 



1 dur + 



\{eff + {p^f- 2pf9^)xf A ■ 

^9t ]dt 



1 + 



1 dMr + 



1 + 



1 dMf^ + 



1 + 



(pi 



ipt + 



Xfdt 



Pt 



^t \ 



)i-W + pt 



Pt 



^t \ 



) Kdt 



then if we set for all t < T 



Pt--- 



Pt 



crt 



and using the bound condition of Q, p,a^,9^, the process is a true martingale. 
Hence we get: 



E 



Qt 



t < T. 



Since Qt = !> we find again the expected result. Moreover L = Z, g^ = —-rxU/^ , 
and gf = 0. Finally, Yt=E and 6 = 0, t < T. 

□ 



Remark 3.8. We have proved that we can find the existence of solution of the triple BSDEs 
using only the VOM. 
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3.3.2 Incomplete market 

In the incomplete market case, the remark 3.8 does not hold true. The VOM depends 
on the dynamics of {QjO'^jO^ , (3). In the particular case where G = F V H"^, we can find 
that the Proposition 3.6 holds true. But in the more general case we 
can not prove the existence of the VOM but we can still characterize the process Q with 
some martingale measure. 

Proposition 3.7. Let consider the incomplete market G = F V IHI"^, then the VOM P 
defines the local martingale measure Q which minimizes the L'^-norm of Z"^ , represents 
the Radon Nikodym density o/P with respect to P on Qt and Zt = K \ZT\Qt\ ■ We find, for 
all t <T, that 



IE miQt] ■ 
Moreover 

Yt = E [^l^\gt] . 

In the more general case, where G = F V H"^ V , we can only prove that there exists a 
martingale measure P such that for all t < T: 

Proof. First step: Consider the case where G = F V and Q be a martingale measure 
for the asset D"^. Let define Z^ its Radon Nikodym density with respect to P on Qt- 
We define the process Z'^ = E Z'^\Qt ■ Using martingale representation Theorem 
1.1, there exists two G-predictable processes and p such that 

dZ^ = Z^_ [ptdM^ + ptdWt] . 

Using Ito's formula, we find: 



d 



^^^P'' l]dMt'^ + {2pt-/3t)dWt + jtdt 



1 + 



3A 



(3.35) 



where jt = {pt — Pt) + . „a + 9t- Since Q is a martingale measure for D we 



get using (2.13) that 

+ pf<^f^f + Pt(^t = 0. 

Hence using this equation, we can find p^ using p and plotting this result on the 
expression of j. We obtain 



AaA\A\'^ ( ,, I fl„ 1 qA„A\A\'^ 

jt = {pt - PtY + 



,2 , (Mi + <ytPt + CJ^O^K) il^t + Ptat + 9fafXi 



(1 + ef){affxf (1 + ef){affx^ + ' 

Let now define 

Pt = Pt-Pt, at = + (1 + ef){aff and h = pt + (Tt(3t + afofX^, 
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then we get: 



3t 



atpt + 2ptbt(Tt + 



at 



at 



p-. + ^) >o, 
at 



(3.35), J > and the fact that the process 
a martingale and ^ G S°°[0,T]), we deduce E 



e 



is a submartingale (since is 



(41 

0T 



> 



since Qt = 1 and 



1. FinaUy, we get for any martingale measure for D"^ that E 



> 



Moreover, if we set pt = then Z is a true martingale measure since 

{e, 6^,6^/3) e 5°°[0,r] X 5°°[0,r] x 5°°[0,r] x BMO and /i,cT^,(T-^ are bounded 
(the process b,a, p and p"^ are bounded). We call P the martingale measure under 
this condition, then E [Z^] = We deduce IP is the martingale measure which 



minimizes the L^-norm of Z and 0+ 



Using the explicit expression of p, we find: 



Pt 



■ — h Pt and pt 

at 



t < T. 



= \- ^t ■ 

at 



Moreover since 



9t 



M^tZt + (1 + 9f) U,^afXt) ^r.^ ^jjA.A 
1- Pt'^t + t^t \ 



at 



+ A + u( 



hot 

at 



Ztpt + u^pfxt 



{i + et)cjfbt 

at 



+ 



then we conclude that Yj = E [-f/jltyt]- Therefore the characterization of the price of 
if) (using Mean- Variance approach) and the VOM in this incomplete case are well 
defined using the vector {S, 9^, 6^,(3) associated to the first BSDE. 



Second step: We consider now the more general case where G = F V V M^. Let 
consider Q be a martingale measure for the asset and let define Z'^ its Radon 
Nikodym density with respect to ¥ onQx- We can define the process 
Using martingale theorem representation 1.1 there exists G-predictable processes p"^, 
p^ and p such that 

dZ^ = Z^_ [pfdM^ + pfdM^ + ptdWt] . 
Using Ito's formula, we find: 



[Zf 




i + ei 



1 dMl + {2pt- f3t)dWt+jtdt 
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where jt = {pt — Pt) + Z^igjA b} i+g'' ^\ 9t- Since Q is a martingale measure for 
we get by (2.13) 

i€{A,B} 

Hence using this equation, we can find using p and p^ and then plotting this result 
on the expression of j. Let first recall a notation: 

«i = ^?+ Yl "'^d bt = Pt + Tt(3t + Yl ^t'^i^t^ 

ie{A,B} iG{A,B} 

SO we obtain: 



B 



i€{A,B} 



+ 



,2 r 



at L 



+ (1 + efXafYxf + 2(pf - - Ai 



+ 26, [(p, - V, + (pf - 0f Vf Af ] . 

Then from the two first terms, we add and remove an additional process to find the 
process a. We get: 



a. 



{pt - I3t) at + -aj + 2ht{pt - Pt)<Tt 
at 



+4(^f)1+(i+w 



+ (1 + 0f )Af 

2{pt-^3t)"^' 



B 



2/ Bn2 {pt 



(1 + 



(1 + 



t J „2 



Finally, we find a more explicit expression of C: 



Ct 



at 



at 



{pf-0f) , btar'' 



+ 



at 



il + 9f)\f{aff[{pt-^) 



Pt 



B _ gB\ 



of 1 + 9P 



It follows that if we set pt - ^t ■= -^^crt and pf - Of := -{1 + 9f)aff^, then we find 
J = and pf-ef = -crf^. Since (|,0^,e^,/3) G S°°[0,T] x5~[0,r] x5~[0,T] x 
BMO and p, a^, and a are bounded then the processes 6, a, p, p"^ and p^ are 
bounded too. Therefore, we deduce that there exists a martingale measure P such 
that 



s<et 



t < T. 



(3.36) 



Moreover we find, for all t < T, that 



gl = ZtPt+ Y 

i&{A,B} 



25 



then 



□ 



Remark 3.9. (About the VOM) 



To identify that P is the VOM in the general case where G = F VH"^ VH^, we should 
prove that j > (as in the first case of the previous Proposition). But from the last 
expression of j, we can not prove that this condition holds true. However, we can 
remark that the assertion of VOM will he justify if one of the following equality is 
satisfied: 



Pt 



1 + 



a 



Apt 



B 



RPt 



1 + 



1 + ef 



or crf{pt - /3t) = CTt 



Pt 



9f 



1 + 



The generalization of the expectation under a a-measure (Yt = K [ip\Gt]) was defined 
by Cerny and Kallsen in [3] p 1512. Moreover, given the solution of the first BSDE: 
(6, e^, e^,P)e 5°°[0, T] X 5~[0, T] X 5~[0, T] x BMO with the constraint O > 6 > 0, 
the martingale Z satisfies (3.36). Since ip is bounded, we conclude: 



\Yf\ < 2E 



Z2 



+ 2E[[V^|2|g;i] <2 



1 

1+ 



Therefore Y G 5°°[0,T] and from representation theorem 1.1, the martingale part M 
ofY given by 

Mt= [ V UidMi + / ZsdWs 
■^0 ie{A,B} 

is BMO. Moreover from Lemma 3.1 in Appendix, since Y € 5°^[0,T], we obtain 
that 9^ and 6^ are bounded. We conclude if the solution of the first BSDE exists 
{e,9^,9^,l3) € 5°°[0,r]x5°°[0,r]x5°°[0,r]xBMO, with the constraint® > 5 > 0, 
that the solution of the second BSDE {Y, U'^,U^,Z) e 5°°[0, T] x5°^[0, T]xS°^[0, T] x 
BMO exists. 



3.4 Proof of Theorem 3.4 

We prove in this part the existence of (9, 6*^, 6*^, /3) in the space 5°°[0, T] x S°°[0, T] x 
5°°[0,T] X BMO with the constraint Q > 6. Moreover, we recall that given the solution of 
this first BSDE, the existence of the second and the third BSDEs is given by Remark 3.5 
and 3.9. 

Note that to prove the existence of {Q,9^,6^ , /3), we do not need the assumption that 
the VOM exists and should satisfied the M2(P) condition (this assumption implies that the 
Radon-Nikodym of the VOM P with respect to P on Gt is non- negative) . Moreover, if 
{e,e^,e^,l3) is defined such that Z is non negative, then it implies that P satisfies the 
i?2(P) condition. 
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In fact, in the general discontinuous filtration, it is difficult to prove that we can find 
{@,e^,9^,/3) solution of the first BSDE such that 9 > (see [23] for more discussions 
about the difficulty). Indeed in the set up of [23], the author makes the hypothesis that 
all coefficients of the asset are F-predictable. This strong hypothesis makes the jump part 
of process equals to zero. In our framework, this hypothesis can not be satisfied since 
the intensities processes are G-adapted. Hence, we deal with splitting method of BSDE 
defined by [18] to prove that, even if the jump part of the process is not equal to zero, we 
can split the jump BSDE in continuous BSDEs such that each BSDE have a solution in a 
good space. The proof is divided in two parts. Firstly, we will give the splitted BSDEs in 
this framework and secondly, we will solve recursively each BSDE. 

First step: Let define, for ah t G [0,r], gt = Qt-gj, 9\ = Gt-Ol for i G {A,B}, Ot = 
^t^l{t<T-4} + ^f^{T-'^<t<TB} and f3t = Q^-f^t. Then, we can define the BSDE (^,0t) 
which is given by: 

det = -ftdt + efdH^ + efdHf^ + ^tdWt with Gt = 1, 

where ft = gt + df^t' + • We also define 

Afc = {{h, ■ ■ ■ h) G (M+)'^ :/!<•••< h}, l<k<2. 

Since we work with the same assumption (density assumption) and notation as in 
[18], then we can decompose Qt and g between each default events: 

@T = 7°l{{0<T<rA}} +7^('^'^)l{T'4<T<rS} +7^('^^i7-'^)l{rS<T} 

and 

ft{et,etJt) = f?{eu9tJt)i{o<t<r^} + fti&t,etJt,r'')i{^A^t<rB} (3.37) 

+ f^{et,etJt,{r^,r''))l{rB<t} (3.38) 

where 7*^ is J^-measurable, is J-'r(8)i3(Afc)-measurable for k = {1, 2}, is V(¥) (g) 
B{R) (g) B(]R)-measurable and f is V{¥) ® B{R) B{R) O B{A''). Moreover, since 
6r = 1 (bounded) (see proposition 3.1 in [18]), we have that the variables ^^{l{^k)) = 1; 
for k = {0,1,2}. 

Let now give the main result of splitting BSDE which is a first step to prove the 
existence of {Q,0^ ,0^ ,^). Let Z(2) = (hjh) G A2 and assume that the following 
BSDE: 

'^0?(Z(2)) = -A'(e?(/(2)),O,42(/(2)),O^it + f (Z(2))dVFi, 0^1(2)) =7\k2)) (3.39) 

admits a solution (0?((/(fc+i))), /3|((^(fc+i)))) € S'^Hh A T,T]) x n'^[l2 A T,T]. And 
that 

deUhk)) = -ft{eHhk))A&t\hkhi)-®t{hk)))^^Hhk))J{k))dt + ^Hhk))dWu 

Qkhk)) = I'ikk)) (3.40) 
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admits ,for k = {0,1}, a solution {e''{l(^k), p^{l^k))) G S°°{[lk AT,T]) x n^[lkAT,T], 
where 1/^) = {h,---lk)- Then (G, P, /3) is given fohowing [18] by: 



©t = ©t l{t<T'4} + ©f (r^)l{^A<j<^i3} + 6t (r^, T'^)l{^s<f}, 
A = /3°1|,<,A| +^,1(t^)1{,a<,<,b} +/32(r^,r^)l{^,<^j, 
07 = e2(r^,t)-ei(r^), 

et = Q\{t)-Q'i{t)- 



(3.41) 



Therefore, to prove the existence of {Q^O"^ ,6^ , j3) we have to prove the existence of 
solutions of the BSDEs (3.39) and (3.40). 

Second step (the recursive approach): We prove recursively the existence of the so- 
lution of these BSDEs. Firstly, we prove the existence of the solution of (3.39) and 
secondly, assuming that the solution of (3.40) exists and satisfies the constraint for 
step /c + 1, we prove the same assertion for the step k. 

1. Let consider the BSDE (3.39): 

0!©?(/(2)) = -/?(©?(/(2)),0,/32(/(2)),/(2))^it + ^j'(/(2))d^^i, ©1(^(2)) = 7'(/(2)). 

Since the coefficient g is given by: 

- 1 2 



5t(©t,^,A) 



©*^? + E.e{A,B}(©< + ^l)KO'A: 



2-i ' 

t 



(3.42) 



from the predictable decomposition of the composition of assets coefficient: 



o"i = o-°l{i<TA} +0-4 (r^)l{^.4<i<^B} +cJt (r^,r^)l{^s<(}, 
= /^°l{t<T'4} + y\{T^)'^{TA<t<TB} + A^?(-^'^,-^^)l{rS<t}, 
'^t ^{t<TA}+crt [T )l{^A<t<^s} + 0-^ (r ,r )l{^fl<j|, 



cr. 



A 



(3.43) 



CJ- 



2,0-1 , 2,1/ Am , 2,2/ A B\-, 

and by our model assumption (see Remark 3.4), we get 



We so find that: 



^t'°l{t<r^} and Af = \^'\T^)l|^^A<^t<^^B^. 



2,1,'_An 



/7(G?(l(2)),0,/32(Z(2)),/(2)) = G?(/(2)) 



/z2(/(2)) ^;9|(Z(2)) 



Using the result of Section 3.3, in the complete market case when 
conclude: 

.2n ^_ ^tili2)) 



E 



t < r,/(2) G A2, 



t£[0,T] 
G = F, we 
(3.44) 
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where the family of processes Z(.) satisfies the SDE given by 



dWt. 



with Zq{1{2)) = 1- Since /i^ and o"^ are bounded, the martingale Mt{l(2)) '■= 
/o i^^^dWs is BMO. We deduce so that Z{1^2)) satisfies the i22(lP) inequality. 
Moreover 7^(^(2)) = 1, we conclude that there exists a constant (5^ > such that 
for all t G [0,r] and 1^2) G A2, Q't{l{2)) > The existence of ^"^(1^2)) 

IS given 

by the martingale part of the process given by (3.44). Moreover since 0^(/(2)) 
is bounded then the coefficient satisfies a quadratic growth with respect to 
/3^(/(2))- Therefore since the terminal condition 7^(/(2)) is bounded, we conclude 
from Kobylanski [19], that P{^2)) is BMO. 

2. Let assume now that there exists a solution which satisfies the constraint for the 
step k + 1. That means that the pair {Q''^'^{l(k+i)), Pt^^{l{k+i))) ^ S°°[k+i,T] x 
BMO and that there exists a non negative constant 5^^^ such that @''^^{l(^k+i) ) — 
(5*^+^. Let us now proved the existence of the pair (©'''(^(fc))) P'^i^k))) ^ <S°°[h, T] x 
BMO at step k: 

d@Ukk)) = -ft (e*=(/(fc)),(e^i(Z(,),t) - QUhk))),^t{hk))J{k)) dt + pHhk))dWt, 

0^a(fc))=7'G(fc))- 
From the decomposition of (3.42), we find: 



+ 



Let consider the processes: 



(^(fc))- 



nt 



k+l,k. 



(^(fc))At"^"^'''(^(fc))0t'^^(^(fc),O; 



k+l. 



mt ■.= aHlik))^tihk))- 



andiVi = m+Kt+mt, dt = \eUiik))WHhk))\ Pt = etHhk)M<'t^'''ihk)y ^t^'''ihk)) 



^2 . fc + l,fc/ 



and Dt = dt + pt- We define 



ft'-- = fli\&Hhk))\J^Jik)) 



Nl 
Dt 



+ 



G 



fc+i 



(^(fe);*) - I©* (^(fc))l 



,k+l,k 



ikk))^ 



(3.45) 



29 



where Ni = n'l+m'l+K'l+2ntmt+2Ktnt+2Ktmt. Since the process G^'^^(/(fc), t) > 
^fc+i y then there exists a non negative constant c > such that pt > c. 
Hence, we obtain: 



e^i(/(fc),t)-|ef(/(fc))| 



A 



k+l,k 



< 



+ 



nt 



+ 



2ntKt 
Pt 



+ |en^w)|A^'''(/(,)) 



+ 



1717 



2ntmt ^ 2mtKt 



dt 



Pt 



+ 



Pt 



Therefore since all processes 0*^"*"^, /i'^, 0-'^+^'''' and X^+^^^ are bounded, there 
exists bounded processes a, b and c such that: 



-ft 



at\Qtihk))\ + btPt{hk)) + 



\Qtii{k))\ 



The coefficient has a quadratic growth and from Kobylanski [19] since the 
terminal condition 7^' is bounded, there exists a pair (B'^(/(fc)), solution of 

the BSDE associated to {f^{l(k))n^Q'(k)))- Moreover, if we consider the BSDE 



dxt = —h(xt, Zt)dt + ZtdWt with terminal condition xt = 7^(^{ 
the coefficient h is given by : 



1)^ 



1 and where 



htixt,Zt) 



-atxt 



Z'2 

^ - btZt, 



tG[0,T], 



we obtain using Proposition 5.11 in [23] that the solution {x,Z) G S°°[0,T] x 
BMO exists. 

Moreover, there exists a non negative constant 5^ such that xt > S^^a.s. Hence 
we conclude using Comparison theorem of quadratic BSDE (see [19]) that > 
—h = h. Finally, the pair of solution (G*^(Z(fc)), associated to (/*^(^(fc)), 7^(^(fc))) 

satisfies 0t(/(fc)) > 5'' > a.s for all t G [0,T],l^k) ^ ^k- Therefore, from 3.45, 
we conclude f'' = It follows that there exists a solution {Q''{l{k)): ^''{^■(k))) ^ 
S°°[lk,T] X BMO associated to if''{l{k)),l''{kk))) such that Q^mk)) > S'' > 
a.s for aU t G [0,r],/(fc) G Afc. 



3.5 Special case and explicit solution of the BSDE 

We conclude by giving an explicit example of our credit risk model which allow us to 
find explicit solution of each BSDEs. We assume G = F V and that the parameters of 
the dynamics of the asset are constant before and after the default time t^. Moreover, we 
assume that the intensity process is given by = A(l — H^). Using theses assumptions, we 
find an explicit solution of the BSDE associated to {g^, Qt) using the splitting approach. 
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Assumption 3.4. The processes /i, o", (T'^,A in (1.4) satisfy the following sassumptions: 

at = (t{H^) = a^l{rA>t} + f^^l{r'4<t}, 



At - A(Fi^) - Al{^A>t}, 

SUc/l t/iat = ((T°)^ + K^A. 

Proposition 3.8. Under Assumption 3.4, there exists a solution of the BSDE associated 
to {g^, G) given by: 



Qt = exp 



-^{T-t) 

K 



1\ 2 



l{r^<i}, t < T. 



Proof. Let first recall that using the sphtting approach developed by [18] , we can write the 
BSDE before and after the default. We obtain 

Gt = G°1|,<,A| + Gi(r^)l{..4<t}, 

1 _ l.Oi , l,li 

9t — 9t ^{t<T-''}+9t i{T^<t}; 

where G" and G^ satisfy the following dynamics: 

dG? 



g\^\el ef,i3^)dt - /3^dWt + xefdt, g" = i, 

= gl'\@lil),0,/3l{l))dt- /3Hl)dWt, QkO = 1 



with 



and gl'\el{l),0,PHl)) 



where / € Ai and Gj (t) — G^ = OfQ^ (see proof of Theorem 3.4 for more details). Using 



Assumption 3.4, setting (3^(1) = 0, we find that gt' (G|, 0, /3/(Z)) 
1, then G^(/) = G^ and we get: 

el = exp 



Since G^(/) 



t < T. 



To find the solution of the first one BSDE, we set (3^ = and from Assumption 3.4 we 
obtain 

We deduce gl'^[Q^, Of, /3^) = — ^t^A. Therefore we find that G*^ satisfies the dynamics: 

l^'dt Q°-l 



Finally, we get G? = exp -^(T - t) 



and we find the expected result. 



□ 



31 



Appendix 



Lemma 3.1. Let consider X and Y two G-predictahle processes such that for i S {A, B}, 
Yr- = X-T^. Then, Xt = Yt on {ri > t) a.s. Moreover, if X^^ < Yr-, then Xt < Yt a.s on 



{n>t). 



Proof. Assume that X and Y are bounded. If X^-. = 1^. , then \Xt — Yt\dHl = and 



= E 



\Xt-Yt\dH, 



E 



\Xt-Yt\Xldt 



Therefore, we have Xt = Yt on (r* > t). Moreover, if X^-. < 5^^, we consider the predictable 
process V defined as Vt = Ytl^Xt<Yt}- Then T/^i = and by using the first part of the 
proof, we obtain Vj = Yj on (r* > t). The general case follows. □ 
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